High-throughput sequencing of black pepper root transcriptome

Background: Black pepper (Piper nigrum L.) is one of the most popular spices in the world. It is used in cooking and the preservation of food and even has medicinal properties. Losses in production from disease are a major limitation in the culture of this crop. The major diseases are root rot and foot rot, which are results of root infection by Fusarium solani and Phytophtora capsici, respectively. Understanding the molecular interaction between the pathogens and the host's root region is important for obtaining resistant cultivars by biotechnological breeding. Genetic and molecular data for this species, though, are limited. In this paper, RNA-Seq technology has been employed, for the first time, to describe the root transcriptome of black pepper. Results: The root transcriptome of black pepper was sequenced by the NGS SOLiD platform and assembled using the multiple-k method. Blast2Go and orthoMCL methods were used to annotate 10338 unigenes. The 4472 predicted proteins showed about 52% homology with the Arabidopsis proteome. Two root proteomes identified 615 proteins, which seem to define the plant's root pattern. Simple-sequence repeats were identified that may be useful in studies of genetic diversity and may have applications in biotechnology and ecology. Conclusions: This dataset of 10338 unigenes is crucially important for the biotechnological breeding of black pepper and the ecogenomics of the Magnoliids, a major group of basal angiosperms.

A comparative transcriptome analysis reveals expression profiles conserved across three Eimeria spp. of domestic fowl and associated with multiple developmental stages

Coccidiosis of the domestic fowl is a worldwide disease caused by seven species of protozoan parasites of the genus Eimeria. The genome of the model species, Eimeria tenella, presents a complexity of 55-60 MB distributed in 14 chromosomes. Relatively few studies have been undertaken to unravel the complexity of the transcriptome of Eimeria parasites. We report here the generation of more than 45,000 open reading frame expressed sequence tag (ORESTES) cDNA reads of E. tenella, Eimeria maxima and Eimeria acervulina, covering several developmental stages: unsporulated oocysts, sporoblastic oocysts, sporulated oocysts, sporozoites and second generation merozoites. All reads were assembled to constitute gene indices and submitted to a comprehensive functional annotation pipeline. In the case of E. tenella, we also incorporated publicly available ESTs to generate an integrated body of information. Orthology analyses have identified genes conserved across different apicomplexan parasites, as well as genes restricted to the genus Eimeria. Digital expression profiles obtained from ORESTES/EST countings, submitted to clustering analyses, revealed a high conservation pattern across the three Eimeria spp. Distance trees showed that unsporulated and sporoblastic oocysts constitute a distinct clade in all species, with sporulated oocysts forming a more external branch. This latter stage also shows a close relationship with sporozoites, whereas first and second generation merozoites are more closely related to each other than to sporozoites. The profiles were unambiguously associated with the distinct developmental stages and strongly correlated with the order of the stages in the parasite life cycle. Finally, we present The Eimeria Transcript Database (http://www.coccidia.icb.usp.br/eimeriatdb), a website that provides open access to all sequencing data, annotation and comparative analysis. We expect this repository to represent a useful resource to the Eimeria scientific community, helping to define potential candidates for the development of new strategies to control coccidiosis of the domestic fowl. (C) 2011 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Transcriptome Analysis of Renal Ischemia/Reperfusion Injury and Its Modulation by Ischemic Pre-Conditioning or Hemin Treatment

Ischemia/reperfusion injury (IRI) is a leading cause of acute renal failure. The definition of the molecular mechanisms involved in renal IRI and counter protection promoted by ischemic pre-conditioning (IPC) or Hemin treatment is an important milestone that needs to be accomplished in this research area. We examined, through an oligonucleotide microarray protocol, the renal differential transcriptome profiles of mice submitted to IRI, IPC and Hemin treatment. After identifying the profiles of differentially expressed genes observed for each comparison, we carried out functional enrichment analysis to reveal transcripts putatively involved in potential relevant biological processes and signaling pathways. The most relevant processes found in these comparisons were stress, apoptosis, cell differentiation, angiogenesis, focal adhesion, ECM-receptor interaction, ion transport, angiogenesis, mitosis and cell cycle, inflammatory response, olfactory transduction and regulation of actin cytoskeleton. In addition, the most important overrepresented pathways were MAPK, ErbB, JAK/STAT, Toll and Nod like receptors, Angiotensin II, Arachidonic acid metabolism, Wnt and coagulation cascade. Also, new insights were gained about the underlying protection mechanisms against renal IRI promoted by IPC and Hemin treatment. Venn diagram analysis allowed us to uncover common and exclusively differentially expressed genes between these two protective maneuvers, underscoring potential common and exclusive biological functions regulated in each case. In summary, IPC exclusively regulated the expression of genes belonging to stress, protein modification and apoptosis, highlighting the role of IPC in controlling exacerbated stress response. Treatment with the Hmox1 inducer Hemin, in turn, exclusively regulated the expression of genes associated with cell differentiation, metabolic pathways, cell cycle, mitosis, development, regulation of actin cytoskeleton and arachidonic acid metabolism, suggesting a pleiotropic effect for Hemin. These findings improve the biological understanding of how the kidney behaves after IRI. They also illustrate some possible underlying molecular mechanisms involved in kidney protection observed with IPC or Hemin treatment maneuvers.

Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq

Schistosoma mansoni is one of the agents of schistosomiasis, a chronic and debilitating disease. Here we, present a transcriptome-wide characterization of adult S. mansoni males by high-throughput RNA-sequencing. We obtained 1,620,432 high-quality ESTs from a directional strand-specific cDNA library, resulting in a 26% higher coverage of genome bases than that of the public ESTs available at NCBI. With a 15 x-deep coverage of transcribed genomic regions, our data were able to (i) confirm for the first time 990 predictions without previous evidence of transcription; (ii) correct gene predictions; (iii) discover 989 and 1196 RNA-seq contigs that map to intergenic and intronic genomic regions, respectively, where no gene had been predicted before. These contigs could represent new protein-coding genes or non-coding RNAs (ncRNAs). Interestingly, we identified 11 novel Micro-exon genes (MEGs). These data reveal new features of the S. mansoni transcriptional landscape and significantly advance our understanding of the parasite transcriptome. (c) 2011 Elsevier Inc. All rights reserved.

Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range

Abstract Background Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain Aw12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. Results To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. Conclusions Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

The use of Open Reading frame ESTs (ORESTES) for analysis of the honey bee transcriptome

Abstract Background The ongoing efforts to sequence the honey bee genome require additional initiatives to define its transcriptome. Towards this end, we employed the Open Reading frame ESTs (ORESTES) strategy to generate profiles for the life cycle of Apis mellifera workers. Results Of the 5,021 ORESTES, 35.2% matched with previously deposited Apis ESTs. The analysis of the remaining sequences defined a set of putative orthologs whose majority had their best-match hits with Anopheles and Drosophila genes. CAP3 assembly of the Apis ORESTES with the already existing 15,500 Apis ESTs generated 3,408 contigs. BLASTX comparison of these contigs with protein sets of organisms representing distinct phylogenetic clades revealed a total of 1,629 contigs that Apis mellifera shares with different taxa. Most (41%) represent genes that are in common to all taxa, another 21% are shared between metazoans (Bilateria), and 16% are shared only within the Insecta clade. A set of 23 putative genes presented a best match with human genes, many of which encode factors related to cell signaling/signal transduction. 1,779 contigs (52%) did not match any known sequence. Applying a correction factor deduced from a parallel analysis performed with Drosophila melanogaster ORESTES, we estimate that approximately half of these no-match ESTs contigs (22%) should represent Apis-specific genes. Conclusions The versatile and cost-efficient ORESTES approach produced minilibraries for honey bee life cycle stages. Such information on central gene regions contributes to genome annotation and also lends itself to cross-transcriptome comparisons to reveal evolutionary trends in insect genomes.

A quantitative view of the transcriptome of Schistosoma mansoni adult-worms using SAGE

Abstract Background Five species of the genus Schistosoma, a parasitic trematode flatworm, are causative agents of Schistosomiasis, a disease that is endemic in a large number of developing countries, affecting millions of patients around the world. By using SAGE (Serial Analysis of Gene Expression) we describe here the first large-scale quantitative analysis of the Schistosoma mansoni transcriptome, one of the most epidemiologically relevant species of this genus. Results After extracting mRNA from pooled male and female adult-worms, a SAGE library was constructed and sequenced, generating 68,238 tags that covered more than 6,000 genes expressed in this developmental stage. An analysis of the ordered tag-list shows the genes of F10 eggshell protein, pol-polyprotein, HSP86, 14-3-3 and a transcript yet to be identified to be the five top most abundant genes in pooled adult worms. Whereas only 8% of the 100 most abundant tags found in adult worms of S. mansoni could not be assigned to transcripts of this parasite, 46.9% of the total ditags could not be mapped, demonstrating that the 3 sequence of most of the rarest transcripts are still to be identified. Mapping of our SAGE tags to S. mansoni genes suggested the occurrence of alternative-polyadenylation in at least 13 gene transcripts. Most of these events seem to shorten the 3 UTR of the mRNAs, which may have consequences over their stability and regulation. Conclusion SAGE revealed the frequency of expression of the majority of the S. mansoni genes. Transcriptome data suggests that alternative polyadenylation is likely to be used in the control of mRNA stability in this organism. When transcriptome was compared with the proteomic data available, we observed a correlation of about 50%, suggesting that both transcriptional and post-transcriptional regulation are important for determining protein abundance in S. mansoni. The generation of SAGE tags from other life-cycle stages should contribute to reveal the dynamics of gene expression in this important parasite.

Modulation of Maternal Immune System During Pregnancy in the Cow

There is a molecular crosstalk between the trophoblast and maternal immune cells of bovine endometrium. The uterine cells are able to secrete cytokine/chemokines to either induce a suppressive environment for establishment of the pregnancy or to recruit immune cells to the endometrium to fight infections. Despite morphological differences between women and cows, mechanisms for immune tolerance during pregnancy seem to be conserved. Mechanisms for uterine immunesuppression in the cow include: reduced expression of major histocompatability proteins by the trophoblast; recruitment of macrophages to the pregnant endometrium; and modulation of immune-related genes in response to the presence of the conceptus. Recently, an eGFP transgenic cloned embryo model developed by our group showed that there is modulation of foetal proteins expressed at the site of syncytium formation, suggesting that foetal cell can regulate not only by the secretion of specific factors such as interferon-tau, but also by regulating their own protein expression to avoid excessive maternal recognition by the local immune system. Furthermore, foetal DNA can be detected in the maternal circulation; this may reflect the occurrence of an invasion of trophoblast cells and/or their fragment beyond the uterine basement membrane in the cow. In fact, the newly description of exosome release by the trophoblast cell suggests that could be a new fashion of maternal-foetal communication at the placental barrier. Additionally, recent global transcriptome studies on bovine endometrium suggested that the immune system is aware, from an immunological point of view, of the presence of the foetus in the cow during early pregnancy.

Identification of sense and antisense transcripts regulated by drought in sugarcane

Sugarcane is an important sugar and energy crop that can be used efficiently for biofuels production. The development of sugarcane cultivars tolerant to drought could allow for the expansion of plantations to sub-prime regions. Knowledge on the mechanisms underlying drought responses and its relationship with carbon partition would greatly help to define routes to increase yield. In this work we studied sugarcane responses to drought using a custom designed oligonucleotide array with 21,901 different probes. The oligoarrays were designed to contain probes that detect transcription in both sense and antisense orientation. We validated the results obtained using quantitative real-time PCR (qPCR). A total of 987 genes were differentially expressed in at least one sample of sugarcane plants submitted to drought for 24, 72 and 120 h. Among them, 928 were sense transcripts and 59 were antisense transcripts. Genes related to Carbohydrate Metabolism, RNA Metabolism and Signal Transduction were selected for gene expression validation by qPCR that indicated a validation percentage of 90 %. From the probes presented on the array, 75 % of the sense probes and 11.9 % of the antisense probes have signal above background and can be classified as expressed sequences. Our custom sugarcane oligonucleotide array provides sensitivity and good coverage of sugarcane transcripts for the identification of a representative proportion of natural antisense transcripts (NATs) and sense-antisense transcript pairs (SATs). The antisense transcriptome showed, in most cases, co-expression with respective sense transcripts.

Diversification of hAT transposase paralogues in the sugarcane genome

Transposons are abundant components of eukaryotic genomes, and play important role in genome evolution. The knowledge about these elements should contribute to the understanding of their impact on the host genomes. The hAT transposon superfamily is one of the best characterized superfamilies in diverse organisms, nevertheless, a detailed study of these elements was never carried in sugarcane. To address this question we analyzed 32 cDNAs similar to that of hAT superfamily of transposons previously identified in the sugarcane transcriptome. Our results revealed that these hAT-like transposases cluster in one highly homogeneous and other more heterogeneous lineage. We present evidences that support the hypothesis that the highly homogeneous group is a domesticated transposase while the remainder of the lineages are composed of transposon units. The first is common to grasses, clusters significantly with domesticated transposases from Arabidopsis, rice and sorghum and is expressed in different tissues of two sugarcane cultivars analyzed. In contrast, the more heterogeneous group represents at least two transposon lineages. We recovered five genomic versions of one lineage, characterizing a novel transposon family with conserved DDE motif, named SChAT. These results indicate the presence of at least three distinct lineages of hAT-like transposase paralogues in sugarcane genome, including a novel transposon family described in Saccharum and a domesticated transposase. Taken together, these findings permit to follow the diversification of some hAT transposase paralogues in sugarcane, aggregating knowledge about the co-evolution of transposons and their host genomes.

A Novel Stress-Induced Sugarcane Gene Confers Tolerance to Drought, Salt and Oxidative Stress in Transgenic Tobacco Plants

Background: Drought is a major abiotic stress that affects crop productivity worldwide. Sugarcane can withstand periods of water scarcity during the final stage of culm maturation, during which sucrose accumulation occurs. Meanwhile, prolonged periods of drought can cause severe plant losses. Methodology/Principal Findings: In a previous study, we evaluated the transcriptome of drought-stressed plants to better understand sugarcane responses to drought. Among the up-regulated genes was Scdr1 (sugarcane drought-responsive 1). The aim of the research reported here was to characterize this gene. Scdr1 encodes a putative protein containing 248 amino acids with a large number of proline (19%) and cysteine (13%) residues. Phylogenetic analysis showed that ScDR1is in a clade with homologs from other monocotyledonous plants, separate from those of dicotyledonous plants. The expression of Scdr1 in different varieties of sugarcane plants has not shown a clear association with drought tolerance. Conclusions/Significance: The overexpression of Scdr1 in transgenic tobacco plants increased their tolerance to drought, salinity and oxidative stress, as demonstrated by increased photosynthesis, water content, biomass, germination rate, chlorophyll content and reduced accumulation of ROS. Physiological parameters, such as transpiration rate (E), net photosynthesis (A), stomatal conductance (gs) and internal leaf CO2 concentration, were less affected by abiotic stresses in transgenic Scdr1 plants compared with wild-type plants. Overall, our results indicated that Scdr1 conferred tolerance to multiple abiotic stresses, highlighting the potential of this gene for biotechnological applications.

Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

Background: Sugarcane is an important crop worldwide for sugar production and increasingly, as a renewable energy source. Modern cultivars have polyploid, large complex genomes, with highly unequal contributions from ancestral genomes. Long Terminal Repeat retrotransposons (LTR-RTs) are the single largest components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their contribution to the genome and transcriptome, however a detailed study of LTR-RTs in sugarcane has not been previously carried out. Results: Sixty complete LTR-RT elements were classified into 35 families within four Copia and three Gypsy lineages. Structurally, within lineages elements were similar, between lineages there were large size differences. FISH analysis resulted in the expected pattern of Gypsy/heterochromatin, Copia/euchromatin, but in two lineages there was localized clustering on some chromosomes. Analysis of related ESTs and RT-PCR showed transcriptional variation between tissues and families. Four distinct patterns were observed in sRNA mapping, the most unusual of which was that of Ale1, with very large numbers of 24nt sRNAs in the coding region. The results presented support the conclusion that distinct small RNA-regulated pathways in sugarcane target the lineages of LTR-RT elements. Conclusions: Individual LTR-RT sugarcane families have distinct structures, and transcriptional and regulatory signatures. Our results indicate that in sugarcane individual LTR-RT families have distinct behaviors and can potentially impact the genome in diverse ways. For instance, these transposable elements may affect nearby genes by generating a diverse set of small RNA's that trigger gene silencing mechanisms. There is also some evidence that ancestral genomes contribute significantly different element numbers from particular LTR-RT lineages to the modern sugarcane cultivar genome.

Trichoderma harzianum expressed sequence tags for identification of genes with putative roles in mycoparasitism against Fusarium solani

The plant pathogen Fusarium solani causes a disease root rot of common bean (Phaseolus vulgaris) resulting in great losses of yield in irrigated areas of the Southeast and Midwest regions of Brazil. Species of the genus Trichoderma have been used in the biological control of this pathogen as an alternative to chemical control. To gain new insights into the biocontrol mechanism used by Trichoderma harzianum against the phytopathogenic fungus, Fusarium solani, we performed a transcriptome analysis using expressed sequence tags (ESTs) and quantitative real-time PCR (RT-qPCR) approaches. A cDNA library from T. harzianum mycelium (isolate ALL42) grown on cell walls of F. solani (CWFS) was constructed and analyzed. A total of 2927 high quality sequences were selected from 3845 and 37.7% were identified as unique genes. The Gene Ontology analysis revealed that the majority of the annotated genes are involved in metabolic processes (80.9%), followed by cellular process (73.7%). We tested twenty genes that encode proteins with potential role in biological control. RT-qPCR analysis showed that none of these genes were expressed when T. harzianum was challenged with itself. These genes showed different patterns of expression during in vitro interaction between T. harzianum and F. solani. (C) 2012 Elsevier Inc. All rights reserved.

Dissecting Phaseolus vulgaris Innate Immune System against Colletotrichum lindemuthianum Infection

Background: The genus Colletotrichum is one of the most economically important plant pathogens, causing anthracnose on a wide range of crops including common beans (Phaseolus vulgaris L.). Crop yield can be dramatically decreased depending on the plant cultivar used and the environmental conditions. This study aimed to identify potential genetic components of the bean immune system to provide environmentally friendly control measures against this fungus. Methodology and Principal Findings: As the common bean is not amenable to reverse genetics to explore functionality and its genome is not fully curated, we used putative Arabidopsis orthologs of bean expressed sequence tag (EST) to perform bioinformatic analysis and experimental validation of gene expression to identify common bean genes regulated during the incompatible interaction with C. lindemuthianum. Similar to model pathosystems, Gene Ontology (GO) analysis indicated that hormone biosynthesis and signaling in common beans seem to be modulated by fungus infection. For instance, cytokinin and ethylene responses were up-regulated and jasmonic acid, gibberellin, and abscisic acid responses were down-regulated, indicating that these hormones may play a central role in this pathosystem. Importantly, we have identified putative bean gene orthologs of Arabidopsis genes involved in the plant immune system. Based on experimental validation of gene expression, we propose that hypersensitive reaction as part of effector-triggered immunity may operate, at least in part, by down-regulating genes, such as FLS2-like and MKK5-like, putative orthologs of the Arabidopsis genes involved in pathogen perception and downstream signaling. Conclusions/Significance: We have identified specific bean genes and uncovered metabolic processes and pathways that may be involved in the immune response against pathogens. Our transcriptome database is a rich resource for mining novel defense-related genes, which enabled us to develop a model of the molecular components of the bean innate immune system regulated upon pathogen attack.

Characterization of the procera Tomato Mutant Shows Novel Functions of the SlDELLA Protein in the Control of Flower Morphology, Cell Division and Expansion, and the Auxin-Signaling Pathway during Fruit-Set and Development

procera (pro) is a tall tomato (Solanum lycopersicum) mutant carrying a point mutation in the GRAS region of the gene encoding SlDELLA, a repressor in the gibberellin (GA) signaling pathway. Consistent with the SlDELLA loss of function, pro plants display a GA-constitutive response phenotype, mimicking wild-type plants treated with GA(3). The ovaries from both nonemasculated and emasculated pro flowers had very strong parthenocarpic capacity, associated with enhanced growth of preanthesis ovaries due to more and larger cells. pro parthenocarpy is facultative because seeded fruits were obtained by manual pollination. Most pro pistils had exserted stigmas, thus preventing self-pollination, similar to wild-type pistils treated with GA(3) or auxins. However, Style2.1, a gene responsible for long styles in noncultivated tomato, may not control the enhanced style elongation of pro pistils, because its expression was not higher in pro styles and did not increase upon GA(3) application. Interestingly, a high percentage of pro flowers had meristic alterations, with one additional petal, sepal, stamen, and carpel at each of the four whorls, respectively, thus unveiling a role of SlDELLA in flower organ development. Microarray analysis showed significant changes in the transcriptome of preanthesis pro ovaries compared with the wild type, indicating that the molecular mechanism underlying the parthenocarpic capacity of pro is complex and that it is mainly associated with changes in the expression of genes involved in GA and auxin pathways. Interestingly, it was found that GA activity modulates the expression of cell division and expansion genes and an auxin signaling gene (tomato AUXIN RESPONSE FACTOR7) during fruit-set.